####### MAURCIO TAKASHI COUTINHO WATANABE ######## Exerccio 9

getwd()
setwd("C:/Users/Maurcio/Documents/aula9")


#### Exerccio Programar 1.: QUE FRIO!

celsius <-function(x)
	{
	diferenca= 5*(x)-160
	celsius=diferenca/9
	return(celsius)
	}

#testando
celsius

#### Exerccio Programar 2. : Anlise exploratria simultnea de duas variveis

require(plotrix)

datagraph <- function(x,y)
	{ 
	x11()
	par(mfrow = c(2,3))	
	boxplot(x,y,col=c("gray","green"),names=c("x","y"), main="Boxplot das variveis")
	multhist(list(x,y),beside=TRUE,col=c("blue","red"),ylab="Frequncia",main="Dados de frequncia")                 
	plot(y,x,col=c("orange","purple"),xlab="",ylab="",main="Grfico de disperso")
	qqnorm(x,col="yellow", xlab="", ylab="", main="Grfico Q-Q de x")
	qqline(x, col="pink")
	qqnorm(y,col="black",xlab="", ylab="", main="Grfico Q-Q de y")
	qqline(y,col="brown")
	par(mfrow=c(1,1))
	estatistica=list(cor(x,y),summary(x),summary(y))
	names(estatistica)=c("Correlao entre x e y", "Sumrio de x", "Sumrio de y")
	return(estatistica)
	}

testando
a = c(1,2,3,4,5)
b = c(3,6,7,9,4)
datagraph(a,b)

###### Programar 3: ndices de Diversidade de Espcies
### Questo feita em sala de aula

diversity <- function(x, modelo="shannon")
	{
	abun=apply(x,2,sum,na.rm=TRUE)
	pi=t(x)/abun
	pi[pi==0]<-NA
	pi2=pi^2
	logpi=pi*log(pi)
		if(modelo=="simpson")
		{
		diversity=apply(pi2,1,sum,na.rm=TRUE)
		}
		if(modelo=="shannon")
		{
		diversity=(-1)*apply(logpi,1,sum,na.rm=TRUE)
		}
	return(diversity)
	}


#testando
inventado <- matrix(rep(10,10),nrow=10)
inventado
diversity(inventado, "shannon")
diversity(inventado,"simpson")